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ABSTRACT 

Recent observations of z > 7 Lya emitters (LAEs) have derived a variety of Lya lu- 
minosity functions (LFs) with contradictory results, evolution or non-evolution from 
z < 6, the epoch after reionization. This could be because most of z > 7 LFs com- 
prise photometric candidates and might include some contaminations. We conducted 
the Subaru Telescope Faint Object Camera And Spectrograph narrowband NB980 
(A c ~ 9800A, FWHM ~ 100A) imaging and spectroscopy survey of z = 7-7.1 LAEs 
to compare its " contamination- free" result with z > 7 photometric Lya LFs previ- 
ously derived. We imaged the Subaru Deep Field and the sky around a cluster MS 
1520.1+3002 and found one LAE candidate, but spectroscopy did not reveal Lya 
though deep enough to detect it. We calculated the expected number of LAEs in our 
survey, using five z — 7 and three z = 7.7 Lya LFs from recent surveys. Seven of them 
are consistent with null detection (O.I^q ^-I.II^'q LAEs) within errors including Pois- 
son statistics and cosmic variance, but average values (0.7-1.1 LAEs) predicted from 
one z = 7 and two z = 7.7 LFs among the seven indicate nearly a single detection. 
The remaining one z — 7 LF predicts 3.0l^ LAEs. As to z = 7, the discrepancy 
likely comes from different LAE selection criteria. For z = 7.7, there are two possibil- 
ities; (1) If z — 7.7 LAEs are somehow brighter in Lya luminosity than lower redshift 
LAEs, z — 7.7 LF is observed to be similar to or higher than lower redshift LFs even if 
attenuated by neutral hydrogen. (2) All/most of the z — 7.7 candidates are not LAEs. 
This supports the decline of LF from z ~ 6 to 7.7 and reionization at z ~ 6-7.7. 

Key words: cosmology: observations - galaxies: evolution - galaxies: high-redshift. 



1 INTRODUCTION 

Lya emitters (LAEs) can be a probe of cosmic reioniza- 
tion, since their Lya emission is absorbed or scattered by 
neutral hydrogen if the universe is not completely ion- 
ized, causing Lya luminosity function (LF) to decline as 
the fraction of n eutral hydrogen in interga lactic medium 
(IGM) increases (|Rhoads fc Malhotral l200ll ) . From z = 3 
to 5.7, the Lya LF was observed not to evolve (e.g., 
lOuchi et al.ll200Sl ). Also, many authors constructed sta- 
tistically large and uniform samples of z = 5.7 and 6.6 
LAEs with some candidates spectroscopical l y confirmed 
(IShimasaku et al.ll2006l; IKashikawa et alJl20oS lOuchi et all 
l2008l . l2010l : lNakamura et al.ll201ll L They found that the Lya 
LF significantly declines from z = 5.7 to 6.6, suggesting 
that the universe could be partly neutral at z = 6.6. The 
decline of the LF was further supported with large spectro- 
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sc opic samp l es ne wly o btained by independent observations 
oflHu et al. (|2010|) and IKashikawa et al.1 ( |201lT ). Meanwhile, 
lOta et all (|20Qg| . l201o"al ) found that Lya LF also declines 
from z = 5.7-6.6 to 7. These studies all imply that neutral 
fraction might increase with redshift at z > 6. 

Moreover, iHu et al l (|2010h noticed that average Lya 
equ ivalent width (EW) is slightly smaller at z = 6.6 than 
5.7. IKashikawa et al.1 (120111 ) also found that Lya EW distri- 
butions of z ~ 3-5.7 LAEs are similar, but EWs of z = 6.6 
LAEs are smaller. On the other hand, some Lyman break 
galaxies (LBGs) are known to show strong Lya emission, 



from 
from 



3 to 6 
6 to 7 



(Stark et al. 2011 


) but suddenly drops 


(lOno et al.1 12012]; 


Pentericci et alj l201ll: 



ISchenker et~aHl2012] ). The lower EW and Lya LBG fraction 
at z > 6 could be due to the rapid evolution of neutral frac- 
tion from z ~ 6 to 7. This is consistent with the idea of the 
partly neutral Universe suggested by the decline of the Lya 
LF. 
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Figure 1. Transmission of the Subaru FOCAS NB980 and the 
Suprime-Cam i' and z' bands (solid curves). These curves include 
the CCD quantum efficiencies, the throughput of the telescope 
and instrument optics, and atmospheric transmission. The OH 
skylines are overplotted (dotted line). 



However, some authors found a fair number of z = 7 
and 7.7 LAE candidates and claim tha t the Lya LF does 
not evolve from z = 5.7-6.6 to 7-7.7 llHibon et all l20ld . 
l201ll , |2012| ; iTilvi et all l20ld ; iKrug et al.ll2012l ). Their re- 



sults contradict late reionization at z < 7.7. Conversely, 
IClement et aD (j201ll ) also surveyed z — 7.7 LAEs but did 
not detect any LAEs and support reionization at z < 7.7. 
Interestingly, all the LFs implying no evolution at z > 6 
are based on photometric samples. If they suffer some con- 
taminations, the conclusion would be different. The best 
way to reveal this is to identify all the candidates by spec- 
troscopy. However, this requires very expensive campaigns 
that observe many candidates spread over different sky loca- 
tions with sufficiently long integration to ensure detections 
or non-detections of Lya. Alternatively, simpler but indirect 
method is to conduct one imaging and spectroscopy survey 
and compare its contamination free result with the photo- 
metric LFs. 

Here, we use our imaging and spectroscopy survey of 
z — 7-7.1 LAEs for this purpose. We did this survey in 
2003, using the Subaru Telescope Faint Object Camera And 
Spectrograph (FOCAS. IKashikawa et~al] I2002D and a nar- 
rowband filter, NB980 (A c - 9800A, FWHM - 100A; see 
Figure [TJ as a groundwork f or our subseque nt la rger z = 7 
LAE surveys in 2005-2010 bv llve et al l i|2006l ) and lOta et ail 
(20081. l2010ah . Since no detailed information about z ~ 7 
galaxies was available in 2003, we selected LAE candidates 
with tentative color criteria and conducted spectroscopy of 
them. In this paper, we re-analyze the NB980 data, refine 
the color criteria based on recent knowledge of z > 7 galax- 
ies, select z = 7-7.1 LAE candidates and see if their spectra 
were taken in 2003 to derive a contamination- free result. 

In Section 2, we describe our imaging data. We perform 
selection of LAE candidates in Section 3. In Section 4, we 
explain the result of follow-up spectroscopy. In Section 5, we 
compare our NB980 survey with previously derived z > 7 
Lya photometric LFs. We conclude in Section 6. Through- 
out we use an (fi TO , f^A, h) = (0.3, 0.7, 0.7) cosmology and 2" 
aperture AB magnitudes, unless otherwise specified. 




Figure 2. NB980-SDF (left) and NB980-MS1520 (right) images 
taken with two CCDs of the FOCAS. The north is up and the 
cast to the left. The field of view is a circle of 6' diameter. 



2 IMAGING OBSERVATION AND DATA 

2.1 Broadband and Narrowband Images 

The NB980 survey ta rgeted the Subaru Deep Field (SDF, 
IKashikawa et alJl2004h and the sky region around a galaxy 
cluster MS1520.1+3002 (hereafter, MS1520). For the SDF, 
broadband BVRi'z' and narrowband NB816 (A c = 8150A, 
FWHM = 120A) and NB921 (A c = 9196A, FWHM = 132A) 
images were taken wi th the Subaru Telescope Suprime-Cam 
l|Mivazaki et al.ll2002f l by the SDF project an d NB973 (A c = 
9755A, FWHM = 200A) image bv llve et all (|2006T ). All the 
SDF images were convolved to have a common point spread 
function (PSF) of O'.^S. Limiting magnitudes at 3<r with 
2" diameter apertures are (B, V, R, i' , 2', NB816, NB921, 
NB973) = (28.45, 27.74, 27.80, 27.43, 26.62, 26.63, 26.54, 
25.47). The i' and z' images of MS 1520 were also taken 
by the Suprime-Cam and have PSFs of 0'.'66 and 0'.'83 and 
limiting magnitudes (2" aperture, 3ct) of 27.33 and 25.74. 

2.2 NB980 Imaging and Data Reduction 

We imaged the SDF and MS1520 with NB980 and FOCAS 
on 2003 May 7-8. The seeing was / .'7-l'.'2. The 600 sec ex- 
posures were dithered with a simple pattern to minimize 
loss of the survey area. First, every time each exposure was 
taken, the pointing was shifted by 2" in the RA direction. 
After repeating this 8 (5) times for the SDF (MS1520), the 
pointing was returned to the original position in the RA di- 
rection but moved by 3" in the DEC direction. We repeated 
this procedure 3 times. The total integration times were 4 
(SDF) and 2.7 (MS1520) hours. 

We redu c ed th e NB980 data in the same manners as in 
iTakata et ai1 ((2003). The dithered exposures were registered 
and combined to produce final stacked images of the SDF 
and MS 1520 (hereafter NB980-SDF and NB980-MS1520) as 
shown in Figure [2] This reduced the fringing, but slight 
residu al remaine d in the stacked images. The standard star 
Hz44 (|Okeil990D was imaged during the observations to cali- 
brate the photometric zeropoints. They were NB980 = 30.27 
mag ADIT 1 for both NB980-SDF and NB980-MS1520. The 
PSFs and 2" aperture 3<r limiting magnitudes of NB980- 
SDF and NB980-MS1520 were l'.'25 and l'.'27, and NB980 
= 24.67 and 24.44, respectively. Assuming the minimum 
detectable rest frame Lya EW of W™* = 20A (see Sec- 
tion [3T2]) s these magnitudes correspond to Lya flux limits of 
Fhya ~ 1.2 and 1.4 x 10 -17 erg s _1 cm -2 or Lya luminosity 
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Figure 3. Redshift evolution of z' — NB980 colors of various galaxies. The colors of E (elliptical), Sbc, Scd and Im (irregular) 
galaxies were calculated using [Coleman. Wu &c Weedmanl dl980T > template spectra. The colors of LBGs/LAEs were calculated using 
the spectrum fx oc \@ for several UV continuum slopes and Lya EWs W£ya' (Left) Colors of LBGs (WJ eB * = 0). They show excess 
of z' — NB980 > 2.5 at z ~ 7-7.1 (horizontal line), but the max imum excess at z ~ 7-7.05 does not depend on 0. (Right) Colors of an 
LBG and LAEs with /3 = —3 are plotted as a representative, as lOno e t al. (201fJ) suggest that z ~ 7 LBGs and z = 5.7 and 6.6 LAEs 
tend to have j3 ~ —3. The LAEs clearly show the excess of z' — NB980 > 2.5 at z ~ 7-7.1, which we adopted as the LAE selection 
criterion (3) in Section [372] This criterion can detect both z ~ 7-7.1 LBGs (W[ es ^ < 20A) and LAEs (W™^ > 20A). 




Figure 4. The 20" X 20" i', z' and NB980 images (from left to 
right) of MS1520-A. The north is up, and the east is to the left. 

limits of Lj^ya ~ 6.7 and 8.3 x 10 42 erg s _1 , respectively. The 
effective area imaged with the FOCAS and NB980 was ~ 56 
arcmin 2 (SDF plus MS 1520). The comoving distance along 
the line of sight corresponding to the redshift 7.0 ^ z ^ 7.1 
probed by the NB980 was ~ 29.1/i _1 Mpc. Hence, we sur- 
veyed a comoving volume of ~ 10 4 Mpc 3 . 



3 LAE CANDIDATE SELECTION 
3.1 Photometry 

We first matched astrometry and pixel scales of all the 
Suprime-Cam images to those of NB980-SDF and NB980- 
MS1520. Also, we convolved PSFs of the z'-band images 
to those of NB980-SDF and NB980-MS1520 to calculate 
z' — NB980 color by measuring z' and NB980 magnitudes 
with the same aperture to select LAE candidates in Sec- 
tion [3]2] Then, we perfor med source detection an d photome- 
try with the SExtractor dBertin fc Arnouts|[l996l ) . The pixel 
scale of the NB980 images is 0'.'1038 pixel -1 . We regarded an 
area larger than 5 contiguous pixels with a flux greater than 
2<r as an object. We detected objects in NB980-SDF and 
NB980-MS1520 and did photometry in other wavebands, 
using the double-imaging mode. The 2" diameter aper- 



ture magnitudes were measured with the MAG_APER param- 
eter and total magnitudes the MAG_AUT0. Finally, combining 
the photometry in all the wavebands, we constructed the 
NB980-detected object catalogs for the SDF and MS1520. 

3.2 LAE Criteria and Candidate Selection 

Figure [T] shows that the NB980 band is located at the red 
side of the z band. If the Lyman break of an LBG is 
redshifted into NB980, it results in a significant excess of 
z' — NB980. If the spectrum has Lya emission, the excess 
is more significant. We used this characteristic to isolate 
2 ~ 7 LBGs/LAEs. We examined the expected z' - NB980 
colors of z ~ 7 LBGs/LAEs and derived candidate selec- 
tion criteria. We first created model spectra of LBGs hav- 
ing the power law continua fx oc A* 3 with several different 
slopes P = —3, —2.5, —2, —1.5, —1, and then model spectra 
of LAEs by adding Lya emission with rest frame EWs of 
W£ya = 20, 50, 100, 150 and 300A. We did not assume any 
specific line profile or velocity dispersion. Instead, we sim- 
ply added the total line flux value to the spectra at 1216A. 
Then, we redshifted the spectra to z = 0-8 a nd applied Ly a 
absorption by IGM, using the prescription of Madaul l|l995h . 

Colors of LBGs/LAEs were calculated using the model 
spectra and transmission curves of z and NB980 and plot- 
ted as a function of redshift in Figure [3] For comparison, 
we also calculated the colors of E (elliptical), Sbc, Scd and 
Im (ir regular) galaxies using the lColeman. Wu fc Weedmanl 
(1980) template spectra. As clearly seen, LBGs (W^ya = 0) 
and LAEs (W^yS ^ 20A) are expected to produce significant 
excess in NB980 against z at z ~ 7-7.1. Also, z ~ 1-3 ellip- 
ticals show modest excess due to the 4000A Balmer break. 
We adopted z' — NB980 > 2.5 as a criterion that selects 
z ~ 7-7.1 LBGs/LAEs avoiding ellipticals. We used the fol- 
lowing criteria to select LBG/LAE candidates. 

NB980 3cr limiting magnitude (1) 



4 Kazuaki Ota and Masanori lye 



Table 1. Photometry of MS1520-A. 



RA(J2000) 


DEC(J2000) 


i' z' NB980 


NB980 


15:22:11.8 


+29:50:40.3 


>28.52 >26.94 23.87 


23.42 



NOTE: The columns 3—5 (6) are 2" aperture (total) magnitudes. 
The limits are la. 



Bands blueward of Lya < 2a, and no object is seen 
visually (2) 
z - NB980 > 2.5 (3) 
NB973 < 2cr limiting magnitude, and 
the object is seen visually (only for the SDF) (4) 

The (1) is the object detection limit. The (2) means nonde- 
tections in B, V, R, i', NB816 and NB921 for SDF and in 
i' for MS1520. No object should be seen by eyes in any of 
these bands since the flux blueward of Lya should be ab- 
sorbed by IGM. This eliminates interlopers such as L/M/T 
type stars and low redshift galaxies such as H/3, [OIII], [Oil], 
Ha or [SII] line emitters. The (3) is the NB980 excess as ex- 
plained above. In calculating the color, if the z' magnitude 
was fainter than 1<t, it was replaced by the \a value. The 
(4) is NB973 < 25.9 meaning a detection in NB973 whose 
bandpass targets z = 6.9-7.1 Lya. We have the NB973 im- 
age of the SDF. Because the NB973 image is deeper than 
the NB980 image, any object detected in the NB980 image 
should be detected and visually seen in the NB973 image. 

We applied the criteria (l)-(4) to the NB980-detected 
object catalogs. Then, we visually inspected the NB980 im- 
ages of selected objects to remove obviously spurious ones 
such as tails of saturated pixels from bright stars, halos of 
bright stars, noises of anomalously high fluxes at the edges 
of CCDs and field of view, and residual of fringing. As a re- 
sult, one object in MS1520 remained (hereafter MS1520-A). 
Table 1 and Figure 4 show its photometry and images. 



4 FOLLOW-UP SPECTROSCOPY 

In our previous NB980 survey in 2003, we had selected more 
candidates with more inclusive criteria z' — NB980 > 1.3 
and 2a limiting magnitudes. Then, we had conducted spec- 
troscopy on 2003 June 22-23 with multi-object slits and FO- 
CAS. We confirmed that MS 1520- A selected in Section EH1 
had been also observed in 2003. At that time, we had used 
VPH950 grism (grating of 1095 lines mm -1 and resolution 
of ~ 2500) with 058 order-cut filter (coverage of 580-1000 
nm) and 0'.'8 slits. Integration time was 4 hours comprising 
six 2400 seconds exposures dithered along the slit by ±1". 
The spectra of the standard star Feige 34 |Okl[l990) had 
been also obtained and used for flux calibration. The data 
re duction was per formed in the same standard manners as 
in llve et alj (|2006h . We inspected the sky-subtracted stacked 
spectrum of MS1520-A and could identify neither Lya emis- 
sion, a UV continuum nor any other spectral features. 

4.1 Possibility of the Candidate Being an LAE 

To see if we had enough depth to detect Lya, we compared 
the sky background RMS of the stacked spectrum with Lya 



flux estimated from NB980 magnitude. If we assume z = 
7 and a rest frame Lya EW as low as W£ya = 20A for 
the most severe case, ~ 76% of the NB980 flux comes from 
Lya line, and total NB980 magnitude of MS1520-A, 23.42 
converts to Lya flux, F^ a = 3.7 x 10~ lr erg s _1 cm~ 2 . 
lOta et all (|2010bh estimated that slit loss of the flux in the 
FOC AS spectroscopy of a z = 6.96 LAE IOK-1 jive et all 
120061 ) is - 35% with a 0'.'8 slit under a seeing of 1". We 
also used / .'8 slits for our spectroscopy, and the seeing was 
also ~ 1". If we apply this slit loss, the flux is F^^, = 
2.4 x 10~ lr erg s~ 1 cm -2 . Meanwhile, we used binning of 4 
pixels (equivalent to O'.M, smaller than the seeing of ~ 1") 
in the spatial direction to extract one dimensional spectrum 
of MS1520-A. Calculating the variance in unbinned pixels in 
the dispersion direction at 9750-9850A (NB980 passband) in 
this spectrum, we estimated the sky RMS to be 8.8 x 10~ 19 
erg s _1 cm ~ 2 A -1 . The FWHM of Lya lin e, for example, of 
a z = 6.6 LAE varies from 5.5 to 14. 6A (|Kashikawa et akl 
120061 : iTaniguchi et al.ll2005l ). If we assume the FWHMs of 
z = 7 LAEs are similar, the Lya flux is F^ = (0.48-1.3) x 
10" 17 erg s" 1 cm" 2 . This is 2-5 times fainter than F^ a , 
deep enough to detect Lya. Thus, MS1520-A is unlikely an 
LAE. 



4.2 Possibility of the Candidate Being an LBG 

Another possible origin of MS1520-A is a z ~ 7 LBG. How- 
ever, we could neither see any faint continuum by eye on 
the 2 dimensional spectrum nor in the 1 dimensional one 
after the 4 pixel binning in the spatial direction. To see if 
we had enough depth to detect the continuum, we compared 
the sky RMS of the stacked spectrum with the UV contin- 
uum flux density estimated from NB980 magnitude. If we 
assume z = 7.05 Lyman break (i.e., all the NB980 flux are 
from the UV continuum at the red half of NB980 and zero 
flux otherwise), total NB980 magnitude of MS1520-A, 23.42 
and 35% slit loss converts to a UV continuum flux density, 
/ajjv = 6.4 x 10" 19 erg s" 1 cm" 2 A" 1 . The sky RMS of 
the 1 dimensional spectrum, 8.8 x 10 -19 erg s _1 cm' 2 A -1 , 
is 1.4 times shallower than /™ v , not enough to detect the 
continuum. Our spectroscopy does not reveal if MS1520-A 
is a z ~ 7 LBG. 

Another way is to see if o ur NB980 imaging was deep 
enough to detect z ~ 7 LBGs. lOuchi et al.l (|2009l ) recently 
detected 22 z ~ 7 LBG candidates in the SDF and the 
GOODS-N fields to the UV luminosity Mvv = -21. If we 
assume that all the NB980 flux is from the UV continuum 
and that 2 = 7, our depths, NB980 = 24.67 (SDF) and 
24.44 (MS 15 20), convert toMuv ~ -22.3 and -22.5. To 
these limits, lOuchi et al. 1 mi) did not detect any LBGs. 
Hence, MS1520-A is unlikely az~7 LBG. 



4.3 Other Possibilities 

The other possible origin of MS1520-A is a late-type star, 
a variable/transient object or a noise. In case of a late- type 
star, the spectrum could show no signal if the continuum is 
fainter than our spectroscopy limit. Because we have only 
one image blueward of Lya (i'-band), the null detection cri- 
terion only on this band might not be strict enough to re- 
move a faint late- type star. In case of a variable/transient 
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object, it could have been fainter than our detection limit 
at the time of i'-band imaging and the spectroscopy while 
it might have been bright at the time of NB980 imaging. 

Meanwhile, MS1520-A is unlikely a low- 2 line emitter. If 
it was, its line flux would be as bright as the FjLa estimated 
in Section |4. II and detectable by the spectroscopy. 

The purpose of this study is to obtain a contamination 
free result for z ~ 7 LAEs, including null detection. We can 
safely conclude that we detect no 2 ~ 7 LAE in our survey 
volume and to our detection limit. We used this information 
to assess z > 7 photometric Lya LFs in the literatures. 



5 DISCUSSION 

We compared our null detection with the expected detection 
number of LAEs estimated from a variety of z — 7 and 7.7 
photometric Lya LFs to check their consistency with the 
current contamination free result. 



5.1 Comparison with z — 7 Lya LFs 

In Tabled we calculated and listed expected detection num- 
ber of LAEs in our NB980 survey volume by integrating or 
interpolating Lya LFs from four recent z = 7 L AE surveys 
(|Otaet alJl200a . l2010al : lHibon et al.ll201ll , l2012l ) to our sur- 
vey limits. The errors include Poisson error for small num- 
ber statistics and cosmic variance (see the footnote of Table 
2 for details). For comparison, we also calculated the ex- 
pected number in the case of no neutral hydrogen attenua- 
tion of Lya emission (i.e., Lya transmission to neutral hy- 
drogen Tf™ 1 = 1) by integrating such a z = 7 Lya LF pre- 



dicted by a recent LAE evolution model of lKobavashi et all 
(|2007h to our sur vey limits. The expected n umbers estimated 
using LFs from lOta et al.l (|200Sl , I2010al 1 and one of the 
iHibon et al.l \20m survey fields (D33) are consistent with 
the null dete c tion. Conversely, the number estimated using 
IHibon et all (120111 ) LF indicates the detection of 3.0+% ;jj 
LAEs. Also, though consistent with null detection within 



the error, the number, 1.1 



+2.2 



estimated from another LF 
(D41) of lHibon et ail ||2012D is larger than those estimated 
from other z = 7 LFs. This implies that LAE candidates 
detected by IHibon et al. I (|201ll . |2012| ) might include some 
degree of c ontamination s . We compa red lOta et al] (2008, 
l2010al ) and IHibon et all |201ll . |2012| ) surveys and investi- 
gat ed what cau s es the discrep ancy. 

lOta et all (|2008l . l2010al ) obtained t he narrowband 
NB97 3 images of the SDF and the SXDS (|Furusawa et all 
2008) fields with one pointing of the Subaru Suprime-Cam 
each, reaching limiting magnit udes NB973 = 24.9 and 25.4 
(5a, 2" aperture). Meanwhile, IHibon et al.l (|2012l ) imaged 
two fields called D33 and D41 with the same NB973 fil- 
ter and one Suprime-Cam pointing each and reached lim- 
iting magnitu des NB973 = 24 3 and 24.7 (5a, 2" aper- 
ture). ThoughlHibon et al l (|2012T i images are sha lower than 
lOta et al. I (|2008l , l2010af T they detected 7 and 7 LAE can- 
didate s in the D33 and the D41, while lOta et all (|2008l . 
l2010al ) detected 1 a nd 3 in the SDF a nd the SXDS (see 
Table©. Moreover. IHibon et all (|201ll ) imaged the COS- 
MOS field with their narrowband NB9680 (A c = 9680A, 
FWHM~ 90A) to detect 6 LAE candidates, though their 



limit and volume are s hallower and smaller than those of 
lOta et all l|200Sl . l2010al). 

Meanwhile, iKobavashi et al.l (|2007l ) z — 7 LF with 
1 predicts the detection of 0.6j| n'£ LAEs even if the 



T, 



IGM 



Lya 



neutr al fraction at z = 7 is 0%. Because IHibon et alj (|201ll . 
120121) LFs predict l.l+"-3.0±§;o detections whether or not 
neutral fractions is 0%, the di ff erenc e in th e exp ected LAE 
numb e r betw een IHibon et all (|201ll . |2012| ) and lOta et all 
l|2008l . l2010al ) cannot be explained by field-to-field varia- 
tion of the degree of Lya attenuation by neutral hydrogen . 
In add ition, since IHibon et al.l (|2012T ) and lOta et"atl l|2008l . 
2010a) surveyed similar and very large volumes, cosmic vari- 
ance and Poisson error are not the cause of difference, ei- 
ther. Moreover, b ecause Ota et al. J2010al) probed to much 
deeper limit than IHibon et aL I l|201ll . l2012l ). if difference in 
dust extinction of all the z — 7 LAEs in different sky fields is 
E(B — V) < 0.06 (equivalent to difference in survey depth), 
dust extinction is also unlikely the reason. For example, 
lOno et al l l|2012D constrained dust extinction of a spectro- 
scopically confirmed z = 7.213 LAE to be E(B - V) ~ 0.05. 
Though it is one example, dust extinction of individual z ~ 7 
LAEs seem to be modest. Thus difference in dust extinction 
among 2 ~ 7 LAEs could be even smaller. One remaininj 



factor that causes the di s crepancy be tween lota et all (|2008l . 
l2010al ) and lHibon et all l|201ll . l2012l ) is different LAE selec- 
tion criteria. 



5.1.1 Effect of Selection Criteria on Detection Number 

IHibon et al.l l|201ll . I2012T I did not consider the redshift evo- 
lution of nar rowband excess of galaxies, while lOta et al.l 
420081 l2010al ) did. The NB9680 and NB973 they used have 
FWHMs of 90A and 200A, and located at the red edge of 
z'-band. This makes z < 7 LAEs/LBGs also selected as can- 
didates. In the case of NB973, z > 6.4 Lya emission with 
EWs of even < 20A can produce si gnificant excess o f 

z' - NB973 > 0.7, according to Figure 3 of lOta et all (|2008l ) 
that plots redshift evolutio n of z' — NB973 colors of model 
LAEs/LBGs generated with lBruzual fc Chariot] (|2003l ) pop- 
ulation synthesis models. The same Figure also shows that 
a 2 = 7 LAE/LBG is expected to have z - NB973 > 1.9. 
We also calculated 2' — NB973 color versus redshift simi- 
lar to Figure 3 in the present paper, using spectra fx oc A^ 
with P = -3 to and = to 300A and found that 

a 2 = 7 LAE/LBG is expected to have 2' - NB973 > 2.2. 

icopically confirmee 
(IOK- l,|Iye_. 1 al. 200(J) has a color of 2' - NB973 > 2.44. 



For exa mple, a spectro scopically confirmed 2 = 6.96 LAE 

~pJT 

lOta et al.l (2008, 2010a) adopted slightly less strict criterion 
2' — NB973 > 1-1.72, considering the effects of photometric 
errors and possible diversity of 2 — 7 LAEs on the color. This 
criterion i n prin c iple se lects z > 6.5 LAEs/LBGs. Hence, 
lOta et al. I (|2008l . l2010al ) imposed an additional criterion, 
a null detection in the narrowband NB921 (A c = 9196A, 
FWHM = 132A) to avoid 2 = 6.5-6.6 LAEs/LBGs. Hence, 
their criteria sel ect 2 > 6.7 LAEs/L BGs. 

Meanwhile, IHibon et all (|2012T l used even more inclu- 
sive criterion 2' — NB973 > 0.65, which selects 2 > 6.4 
LAEs/LBGs, and they did not have NB921 images to avoid 
2 = 6.5-6.6 LAEs. Their sample include 4 and 6 candi- 
dates with 2' - NB973 ~ 1-1.5 in the D33 and the D41 
fields. Some of t hem m ight be 2 = 6.5-6.6 LA Es, which 
lOta et all (|2008l . l2010al ) did not select. In fact, lOta et all 
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Table 2. Expected detection number of LAEs in the NB980 survey estimated from z = 7 and 7.7 Lya LFs. 



Authors 


Volume" 


Lya flux limit" #LAE candidates'* 


Expected #LAEs in 


NB980 survey 6 




(10 4 Mpc 3 ) 


(10 — 17 erg s — 1 cm~ 2 ) authors detected 


SDF 


MS1520 


SDF+MS1520 


2 = 7 LAEs 


This study 


1.0 


1.2-1.4 











Ota et al. (20081 


32 


1.5* 


1 


3+ 19 


2+ 1 ' 8 


3+ 18 


Ota et al. f 2010a) 


30 


0.97 6 


3 


u - i -o.i 


0.07l ; 8 7 


o i+i- 8 


Hibon et al. (2011) 
Hibon et al. (20121 D33^ 
Hibon et al. ('2012') D41^ 
Kobavashi et al. ( 20071 


7.2 
40 
43 


1.8 C 
2.7 b 
1.9 b 


6 
7 
7 


2 0+ 2 ' 8 

5+ 1 - 9 

1 1+ 23 

4+1-9 


1 5+ 1 - 9 

1- S -1.5 
2+ 1 - 8 

5+ 1 ' 9 
u -°-0.5 
3+1.8 


3 0+ 32 

5+ 1 ' 9 
U -°-0.5 

1 1+ 2 - 2 

6+ 2 - 7 
u -°-0.5 


z = 7.7 LAEs 


Hibon et al. (20101 
Tilvi et al. (20101 


6.3 
1.4 


0.83 
0.60 


7 
4 


3+1.9 
7+ 2 - 6 

u -'-0.6 


2+ 18 
4+1-9 


3+1.8 


Krusr et al. (20121 


2.8 


0.80 


1 


5+ 19 
u -°-0.5 


4+1-9 





a The survey volume and limit of the papers in the column 1. 'Lya fluxes converted from their limiting NB973 magnitudes (5<r, 2" 
aperture), assuming a rest frame Lya EW of 20A and z = 7. c The Lya flux directly converted from the NB9680 magnitude of their 
faintest LAE candidate. d The number of LAE candidates the authors detected in their survey volumes and to their survey li mits. e The 
expec ted detection number of LAEs in our NB980 survey estimated using the Lya LFs from the papers in the column 1 . For lOta et ah! 
(20081, wc integrated the inferred z = 7 Lya LF i n Figur e 10 of their paper. lOta et al.l l l2010al l Lya LF is based on LLya directly converted 
from Lya fluxes of their can didates. lOta et al l l|2010rj l estimated that ~ 77% of NB973 total flux of a z = 6.96 LAE (IOK-1) is from 
Lya. Hence, we interpolated Ota et al. ll2010all Lya LF attenuated by 0.77 X L^ v a to our NB980 survey limits to obtain the expected 
nu mber of LAEs. Forlffibon et al.l (|2010l . T201ll . l2012f ). we integrated the lSchechten l|l976r i LFs best-fitted to their obse rved data by them . 
For lKobavashi et al.1 d2007t) . we integrated the T 1 ^ 1 = 1 Lya LF predicted by their model. For lTilvi et~al] | |2010]1 and lKrug et alj <2012l) , 
we interpolated their Lya LFs to our NB980 survey limits. The errors include Poisson errors for small number statistics dGehrelsHl986T) 
and cosmic variance a v . We combined these errors quadratical ly and if the lower limit was a negative value, we adjusted it to zero. We 
calculated a v 's using the bias b = 3.4 obtained for LAE s by lOuchi et al.l ||2005| 1 and the dark matter halo variances ctdm's at z = 6 
predicted by the analytic model of lSomerville et~al 1 (|2004l l and our survey volumes. The a v 's are ~ 48% (~ 54%) for the survey volumes 
equivalent to two (one) FOCAS field of views, f The Schechter parameters <E>* and L* for LFs derived from the D33 and the D41 fields 
in Table 4 in I Hibon et al] l|2012l l paper was found to be mistakenly listed. They are switched with each other (P. Hibon 2012, private 
communication). We used their Schechter LFs with correct parameters to estimate the expected LAE numbers. 



(2008) reported that if they did not impose null detection in 
NB921, they selected a spectroscopically confirmed z — 6.6 
LAE IOK-3 with z ' - N B973 ~ 1.41 as a candidate. Be- 
sides, iHibon et ai] (|201lh adopted NB9680 - z < -0.75 
without examining redshift evolution of this co lor of LAEs. 
For the same reason as the case of NB973, IHibon et al] 
|201ll l sample might potentially include s ome lower redshift 
LAEs/LBGs. If IHibon et all (|201ll . l2012h candidates really 
include some such contaminations, their real z — 7 Lya LFs 
after corrected for the contaminations could show evolution 
from z ~ 6 and predict more consistent number of LAEs in 
our NB980 survey. 



5.2 Comparison with z — 7.7 Lya LFs 

In Table [2l we also listed the expected numbe r of LAEs in 
our survey estimated using 3 z = 7.7 Ly a LFs (|Hibon et al] 
l20ld ; lTilvi et al]|2010l ; lKrug et al]|2012l ) under the assump- 
tion that LF does not evolve from z = 7 to 7.7. All the ex- 
pected numbers are consistent with null detection within er- 
rors. Nonetheless, the average value 0.7-0.8 of the expected 
number [0.7^- O.8t nj LAEs) estimat ed from the LFs of 
iTilvi et al.1 |2010l l and lKrug et al.1 (|2012T ) is close to a single 
detection that means evolution of LF from z = 7 to 7.7. 
Actually, these authors and even IHibon et al] (|2010h claim 
that if their candidates are all real z — 7.7 LAEs, their LFs 



imply a trend opposite to the decline of Lya LF from z = 5.7 
to 6.6 and 7 found by Kashikawa et al. l201ll),IOu~chi et al] 
(l2010h and lOta et al.1 (|200Sl , l2010al ). ITilvi et al.l (|20ld ) and 
iKrug et al l (|2012l ) also performed Monte Carlo simulations 
to estimate the expected number of LAEs in their surveys 
and concluded that even if 1 or 2 candidates in their sam- 
ples are real z — 7.7 LAEs, it implies no evolution of Lya 
LF from z — 6.5 to 7.7. There are two possible scenarios 
that explain the situation. 

One is that z = 7.7 LAEs might be somehow brighter 
than z — 6.5 LAEs in Lya luminosity. In this case, z = 7.7 
LF is observed to be similar to or higher than z = 6.5 
LF even if attenuated by neutral hydrogen. There are two 
factors potentially supporting this. First, stellar population 
studies of z ~ 6 massive (~ 10 10 M Q ) LBGs suggest that 
they could have had substantially high star formation rates 
in the past at z ~ 7-8 to ass emble their high masses and 
could be very bright at z > 7 i|Yan et al]|2006l ; lEvles et ai] 
2007). Also, narrowband selected LAEs tend to be the lower 
mass and younger age extension of LBG population and 
might be at the earlier stage of t he evolution of LBGs 
l|Lai et al] 120081 ; iPirzkal et al]|2007l ). If z = 7.7 LAEs are 
the progenitors of z ~ 6 mas sive LBGs, they m ight be in- 
trinsically very bright. Second. lrlaves et al] l|201ll ) estimated 
Lya escape fractions faZ a at 2 ~ 0-7.7 using observed data 
in literatures and found that increases with redshift 
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when they fitted a power law to feZ a 's at z ~ 0-6. If /^c" 
is higher in LAEs at z — 7.7 than 6.5, z = 7.7 LAEs are 
observed to be brighter. 

Anoth e r scen ario is that all t he LA E candidates of 
Tilvi et all j2010h and iKrug et all (|2012h and all/some of 
Hibon et all l|201Ch LAE candidates are not z — 7.7 LAEs. 
This agree s with the decline of Lya LF from z = 5.7 to 6.6- 
7 found b- 
lOta et al.1 
other z = 

not detect any candidates. As mentioned earlier, the dror. 



Kashikawa et al.1 (|201lj V [Ouchi et alj (|2010h and 

20081 . 12010 j )7 io is supporte d by an- 

7.7 LAE survey bv lClement et al.1 l|201ll ) who did 

2P 

in fraction of Lyq em i tting LBGs at z > 6 |Ono et al .112013 : 
iPentericci et al. I l201ll : ISchenker etal]|2012h could be also a 
support. 



6 CONCLUSION 

We conducted the Subaru FOCAS NB980 imaging and spec- 
troscopy survey of z = 7-7.1 LAEs in the SDF and MS1520 
and detected no LAEs to a 3a Lya flux limit of ~ 1.4 x 10 -17 
erg s _1 cm" 2 in a comoving volume of ~ 10 4 Mpc 3 . 

We estimated the expected number of LAEs in our sur- 
vey from 5 z — 7 and 3 z — 7.7 photometric Lya LFs. 
Seven agree with null detection within errors but average 
LAE numbers predicted by one z = 7 and two z — 7.7 LFs 
among the seven indicate nearly a single detection. The re- 
maining one z = 7 LF predicts a detection of 3.0+2 LAEs. 
For z — 7, the discrepancy likely comes from different LAE 
selection criteria. For z = 7.7, there are 2 possibilities. (1) If 
z — 7.7 LAEs are brighter than lower- z LAEs, z — 7.7 LF is 
similar to or higher than lower- z LFs even if attenuated by 
neutral hydrogen. (2) If all the z — 7.7 candidates are not 
real, Lya LF declines from z ~ 6 to 7.7. 

Evaluation of z > 7 Lya LFs in this study relies on 
an indirect method. The direct and best way is secure spec- 
troscopy of all the candidates. Future powerful telescopes 
such as JWST or TMT will facilitate this and reveal the 
real nature of Lya LFs at the epoch of reionization. 
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